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ON THE DISCRETIZATION OF BACKWARD DOUBLY STOCHASTIC 

DIFFERENTIAL EQUATIONS 

OMAR ABOURA 

Abstract. In this paper, we are dealing with the approximation of the process (Xt,Yt, Zt) 
solution to the backward doubly stochastic differential equation (BDSDE) 



X s 



x+ f b(X r )dr+ f a(X r )dW r , 
Jo Jo 

4>(X T )+f f (r,X r ,Y r ,Z r )dr + ( g(r,Xr,Y r ,Zr)dBr- [ 

J s J s J s 



After proving the L 2 -regularity of Z. we use the Euler scheme to discretize X and the Zhang 
On ■ approach in order to give a discretization scheme of the process (Y, Z). 

Oh 

^ . 1. Introduction 

Since the pioneering work of E. Pardoux and S. Peng |PP92j . backward stochastic differential 
equations (BSDEs) have been intensively studied during the two last decades. Indeed, this notion 
has been a very useful tool to study problems in many areas, such as mathematical finance, 
stochastic control, partial differential equations; see e.g. |MY99| where many applications are 

■ described. Discretization schemes for BSDEs have been introduced and studied by several authors. 
[ The first papers on this topic are that of V. Bally |Ba97] and D.Chevance |Ch97j . In his thesis, 

Zhang made an interesting contribution which was the starting point of intense study among which 
the works of B. Bouchard and N.Touzi |BT04j . E.Gobet, J.P. Lemor and X. Warm[GLW05J,... The 
notion of BSDE has been generalized by E. Pardoux and S. Peng [PP94| to that of Backward 
Doubly Stochastic Differential Equation (BDSDE) as follows. Let (O, J 7 , P) be a probability space, 

■ T denote some fixed terminal time which will be used throughout the paper, (Wt) 0<t<T and 
(-^t)o<t<T ^ e two independent standard Brownian motions defined on (17,^-", P) and with values 
in R d , and R respectively On this space we will deal with two families of a- algebras: 

T t :=J^ t W^ T WM t ^:-^ >^tV^ H t = J% T V T* T V jV, (1.1) 

where Tf T := a (B r - B t ;t <r < T), T™ t := a (W r - W ; < r < t) and TV denotes the class of P 

null sets. We remark that {iFtj is a filtration, (Ht) is a decreasing family of u-albegras, while (J-j) is 
neither increasing nor decreasing. Given an initial condition x £ R d , let (Xt) be the d-dimensional 
diffusion process defined by 

X t = x+ [ b(X r )dr+ f a{X r )dW r . (1.2) 
Jo Jo 

Let £ G L 2 (VL) be an R d -valued, JFr-nreasurable random variable, / and g be regular enough 
coefficients; consider the BDSDE defined as follows: 

Y s = £ + J f(r,X r: Y r ,Z r )dr 

+ f g(r,X r ,Y r ,Z r )dB r ~ f Z r dW r . (1.3) 

J s J s 
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In this equation, dW is the forward integral and dB is the backward integral (we send the reader 
to |NP88j for more details on backward integration). A solution to (|1.3p is a pair of real- valued 
process (Y t ,Z t ), such that X t and Y t are (Ft) for every t G [0,T], such that (|1.3|) is satisfied and 

»T 

E( sup |F s | 2 )+e/ |Z s | 2 ds<+oo. (1.4) 

In [PP94 Pardoux and Peng have proved that under some Lipschitz property on / and g which 
will be stated more precisely in section [21 (|1.3[) has a unique solution (Y, Z). 

The aim of this paper is to study the discretization of a Backward Doubly Stochastic Differential 
Equation For the sake of simplicity, as in Zhang's paper |Z04j . we assume that Y and Z are 
real-valued processes. The extension to higher dimension is cumbersome and without theoretical 
problems. This discretization scheme of (Y, Z) is motivated by the link between (|1.3|) and the 
following backward stochastic partial differential equation when £ = 4>(Xt) for a regular function 



u(t,x) = 4>(x) + / ^£it(s, x) + f (s, x, u(s, x), V«(s, x)a(x)) j 



ds 

g (s, x, u(s, x), Vu(s, x)a(x)) dB s , (1-5) 
where C is the differential operator defined by: 

1 d d 2 d d 

Cu(t,x) = - ]T (aa*) itj (x) gx gx u(t,x) +^2bi(x)—u(t,x). 

— l 1 3 i— 1 

The paper is organized as follows: first we prove the L 2 -regularity of Z in section O This a crucial 
step in order to the scheme using Zhang's method, which is done in section [31 Finally, a numerical 
scheme is described in the last section. To ease notations, we set O r := (X r , Y r , Z r ) for r G [0, T]. 
As usual, we denote by C p a constant which depends on some parameter p, and which can change 
from on line to the next one. Finally, for some function h(t, x,y,z) defined on [0, T] x R d x R x R, we 
let d y h(t, x, y, z) (resp. d z h(t, x, y, z)) the partial derivatives of h with respect to the real variable 
y (resp. z), while d m h(t, x, y, z) will denote the vector (d Xi h(t, x, y, z), i = 1, • • • , d). 



2. Regularity properties 

In this section we give some regularity properties of the process X, Y and Z . 

The following assumptions which ensure existence and uniqueness of the solution will be in 
force throughout the paper. For every integer n > 1, let AI 2 ([0, T],R") denote the set of R n - 
valued jointly measurable processes (cpt, t £ [0, T]) such that tpt is ^-measurable for almost every 
t and E J Q \(p t \ 2 dt < +oo. 

Assumption 1 (for the forward process X). The maps b : R d -> R d and a : R d -> R dxd are of 
class C^. 

Assumption 2 (for the backward process (Y,Z)). Let f : Q x [0, T] x R d x R x R -> R, g : 

O x [0, T] x R d x R x R — > R be such that f and g are jointly measurable, for every (x, y, z) € R d+2 , 
f(.,x,y,z) and g(.,x,y, z) belong to M 2 ([0,T],R), and such that: 

(i) There exist some nonnegative constants Lf,L g and a constant a G [0, 1) such that for every 
cj efl,t,t' e [0, T], x, x' G R d , y,y' el and z,z' G R 

|/ (t, x, y,z)-f (f, x', y\ z')\ 2 < L f (\t - t'\ + \ x -x'\ 2 + \y- y'f + \z- z'f) , 
\g (t,x,y, z) -g(t',x',y', z')\ 2 < L g (\t - t'\ + \x - x'\ 2 + \y- y'f^j +a\z- z'\ 2 , 

(ii) For all s G [0, T] f (s, .) and g(s, .) are of class C 3 with bounded partial derivatives up to 
order 3, uniformly in time. 
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of class C 2 and that for 



(iii) For a function h(t,x,y,z), set h(t, 0) := h(t, 0,0,0). Then 

sup |/(r, 0)|+ sup \g(r, 0)| < oo. 

re[0,T] re[0,T] 

Assumption 3. Suppose that £ := c/)(Xt) for some function <p : M. d - 
every w£fl, 

sup |<9 2 5 z)\ < 1. 

2.1. Some classical properties of the forward process X. We at first recall without proof the 
following well known results on diffusion processes. Define the R dxd - valued process (VX t ) 0<KT 
by: 



VX 



(— 

\dx. 



XI, i,j = !,••• 



Then VX t is an invertible d x d matrix, solution to a linear stochastic differential equation with 
coefficients depending on X t . Furthermore, the assumptions on the coefficients a and b yield the 
following classical result: 

Proposition 2.1. (i) For all p > 1, there exist a constant C p > such that for all t,s £ [0,T]: 



IX - X 



,2p 



(vx t ) '-(vi.r 1 



2p 



<C p |t-*|: 



(zij For o/Z p £ [l,+oo[, </iere errisi a constant C p > smc/i i/iai 



E( sup |X| 2 f+ sup (VX)" 1 

V i£[0,T] t€[0,T] 



2.2. Time increments of Y and F 2 -regularity of Z. The following lemma provides upper 
bounds for time increments of Y. 



Lemma 2.2. Set £ = <P(Xt) for some function 4> '■ R d — > R &e o/ class C\. Then we have 

(i) For all p > 2, there exist a constant C p > depending on T such that for all t,sG [0, T] 

E\Y t -Y s \ p <C p \t-s\i. 

(ii) For allp>\, there exist a constant C > such that 

sup E\Z r \ 2p < C. 

0<r<T 



(2.1) 
(2.2) 



Notice that the inequality (|2.ip is different from equation (2.11) in |Z04j . 

Proof. We at first prove (ii). Let (VYt) 0<t<T = {d x Y t )o<t<T denote the real-valued process defined 
by differentiation of Y as function of the initial condition x of the diffusion process (Xt). We recall 
the following representation of Z (see |PP94j Proposition 2.3): 

z t = vr t (vx)"V(x). (2.3) 

where (VF t , VZ t ) satisfies the linear BDSDE with the forward process (X t , VX t ) and the evolution 
equation: 

VY t =^'(X T )VX T + / (f x (r, e r )VX r + f y (r, G r )VY r + f z (r, Q r )VZ^dr 

(2.4) 
(2.5) 



+ J [g x (r,G r )VX r + g y (r,e r )VY r + g z (r,e r )VZ r )dB r - \ VZ r dW r 
By E.Pardoux and S.Peng [PP94] page 217, we deduce 



E sup |Vy t |" ) < oo. 

■ 0<t<T 



Then Holder's inequality and Proposition 12.11 yield 



E|Z t | 2p < (ElVF, 



1 6 p 



E 



(VX)" 



6p 
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This concludes the proof of (ii) . 

(i) Suppose that s < t, then using (|1.3p . we deduce that 



\Y t -Y s \ p < C v 



f(r,e r )-f{r,X r ,Y r ,0)dr 



C 



-Cm 



g(r, <~> r )dB r 



f(r,X r ,Y r ,0)dr 



Recall that Tt and Jit have been defined in (jl.ip . The process (J ' Z r dW r , < t < T) is a {^Ft)~ 

martingale, while the process (J^ g{r, ® r )dB r , < t < T) is a backward martingale for (7it). 
Hence, the Burkholder-Davies-Gundy and Holder inequalities yield 



E\Y t ~Y s \ p < C^t-s^El \f(r,X r ,Y r ,0)\ p dr + C p E 



\Z r \ dr 



+C P E\^J \g(r,Q r )\ dr) + C P E [ \ \Z r \Ur 
Assumption [2] (i) and (ii), Proposition 12. II and (|1.4|) yield 

E 



(2.6) 



\g(r,9 r )rdr) <C P J£(J \g(r,Q r )-g(r,0)\ dr) + C P EU \g(r,0)\* dr 

t ~ t 

<C p \t-s\S + C p eN (\X r \ 2 + \Y r \ 2 ) d r y +C p eN \Z r \ 2 dr 



<CJt-s\2 + C P 



\Z r \ dr 



(2.7) 



Similarly, 



E \f(r,X r ,Y r ,0)\Pdr<C p E \f(r, 0)\ p dr + C p E / \f (r, X r , Y r , 0) - / (r, 0)| p dr 



<CpE \f (r,0)| p dr + C p E(\X r \ p + \Y r \ p ) dr < C p \t — s\. 

J S J S 

Hence, the inequalities (|2.6p ~ (j2.8p imply 

E\Y t -Y s \ p <C p \t-s\i +C p eH \Z r \ 
Using Holder's inequality and (|2.2j) we conclude the proof of ()2.1 



(2.8) 



' dr 



□ 



Since equation (|2.4[) proves that the pair (VF, VZ) is the solution of a BDSDE with forward 
process (X,VX) G L p for every p G [1, +oo[, we deduce from (|2.ip that for every function <f> : 
R d -> R of class Cl, we have for < s < t < T and p G [1, +oo[: 

E\WY t -VY s \ p <C p \t-s\i, (2.9) 

for some constant C p > 0. We now establish some control of time increments of the process Z, 
following the idea of J.Zhang [Z04] . 

Theorem 2.3 (L 2 -regularity of Z). There exists a non negative constant C such that for every 
subdivision tt = {to = < t\ ■ ■ • < t n = T} with mesh \tt\, one has 

rti 



Ki<n Jt '- 1 



z t -z 



U-i 



\ 2 + \z t - z t j 2 ) dt < cm 



(2.10) 



Proof. Using the representation of Z as a product, we deduce 

Z t - Z ti = VY t (VXt)- 1 a (X t ) - VY U (VA^)" 1 a (X ti ) 
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Then, 



\z t -z t f <3|vr t -vy ti | 2 



3 IVY, 



(vx,)" 1 \<j(X t )\ 2 
{VXtY 1 -{vx u )- l] \\o{x t )\ 2 

(VX,)- 1 2 \a(X t )-a(X ti )\\ 



+ 3|vr t 

To conclude the proof, we use Holder's inequality, Proposition 12. II and (|2.9 

Theorem 12.31 immcdiatly yields the following 
Corollary 2.4. 



□ 



\Z r - Z t .\ 2 dr < CM 



3. The discretization of (X, Y, Z) 

3.1. Discretization of the process X: The Euler scheme. We briefly recall the Euler scheme 
and send the reader to [KP99 for more details. Let 7r := {to = < t\ < ... < t n = T} be a 
subdivision of [0,T]. We define the process X£ , called the Euler scheme, by 



X? = XI + f b (XI) ds+ f a(Xl)dW s 

J to J to 



where := max{ii < s}. The following result is well known: 

Proposition 3.1. There exists a constant C > such that for every subdivision it, 

maxElXt. - X*? < C\ir\, E \X r - X? | 2 dr < C\ir\ 2 . 
i Ju-i 

3.2. Discretization of the process (Y,Z): The step process. In this section, we construct 
an approximation of (Y, Z) using Zhang's approach. 

Let 7r : to = < ... < t n = T be any subdivision on [0, T]. Set Gt = Q\ f° r U-i < £ < ti, where we 
let 

g\ :=o-{W r -W ;0<r<t)Vo- (B r - B^^U-i < r < T) , t*_ x <t<t h 
and define the (<7t)-adapted process (Y^ , Z^) Q<t<T recursively (in a backward manner), as follows: 
Set Y t n n = <j>(X?J, Z^ 1 = 0; for % = n - 1, 0, let 



At 



i+i 



ZZdr 



and for i = n, 1, let 

At i = t i -t i - 1 ,AB ti =B u -B u _„& K t ! 1 := (xi,Y t «,Z?: 1 ) , 



Y t « = Y* + f{t i ,e^) AU + g{ti,Q% 1 ) AB t*-J Z rdW r , Vie (3.1) 
Note that the equation (|3.1|) is not a BDSDE in the sense of [PP94] : however, we have the following: 

Proposition 3.2. For every i = l,...,n, there exists a process (Y* , Z t 7r ) t6 [j._ 1 j.) adapted to the 
filtration (Q t ,U-i < t < U), such that &3. 1\) holds. Furthermore, Y±. £ Tt^. 

Proof. The proof is similar to that in |PP94j page 212 and relies on the martingale representation 
theorem. Fix an integer i > and suppose that the processes (Y*) and (Z£) have been defined for 
t>ti, (£/ t )-adapted, and that is T tk -measurable for k = i, 
the process defined by 



,n. We denote by (M|) te[ti i ti] 



M\ := E (y% +f(ti, ej 1 ) *ti + g (u, O^/ 1 ) AB U 



H—l — t — H 
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By the martingale representation theorem, there exists a (pl,ti-i < t < ii)-adapted and square 
integrable process (^,^-1 <t<u) such that for U-i < t < U, M\ = M l U i + J* N l s dW s . 
Therefore, NP t = M\. - N\dW s . Clearly, Q\. contains T tl , X£ is C T ti measurable and 
is T ti -measurable; hence 

Ml =Y£ + f (U, e^; 1 ) AU + g (u, O*: 1 ) AB U . 



Furthermore, note that G\ i _ 1 = ^U-d so that M*. is -measurable. This completes the proof 
by setting: Y t v = Mj, Z? = Nj for t l _ 1 < t < t % . □ 

Before stating the main theorem of this section, we introduce the following 

Definition 3.3. Let k > 1 be a constant. The subdivision n is said to be K-uniform if kAU > \tt\ 
for every i G {1, n} . 

The main example of a K-uniform subdivision is a uniform subdivision (i.e. for all i, Ati = \tt\) 
where k — 1. The following lemma gives an upper estimate of Z ti — Z^' . 

Lemma 3.4. For any i = 0, ...,n — 1, any K-uniform subdivision tt and (3 > we have: 



Ai,E 



Hi 



< k (1 + 0) j \Z? - Z r \ 2 dr + K (l + /T 1 ) j \Z r - Z t f dr 



Proof. For any i = 0, ...,n — 1, Z ti is Tt t -measurable, and Ati < |tt| < kA^ + i; thus 

2 



ALE 



V 1 



At, 



< 



(At i+ i) 

K 



r3 



AUE 



At l+1 



ZZdr 



At l+1 



-E 



3 

ti+i 



(Z u -Z«)dr 



< kE 



\Z U -Z?fdr. 



where the last step is deduced from Schwarz's inequality. Using the usual estimate \Z ti — Z*\ 2 < 
(1 + @)\Z* - Z r \ 2 + (1 + /3 _1 )|Z r - Z u \ 2 , we conclude the proof. □ 

The following theorem is the main result of this section. It proves that as \tt\ — > 0, (Y^ , Z 71 ") 
converges to (Y, Z). 

Theorem 3.5. Let it be a n-uniform subdivision with sufficiently small mesh a < -, let <f> E C 2 
and £ = 4>{Xt). Then we have 



max E \Y U - Y?? + E / \Z r - Z?\ 2 dr < C\ir\. 

0<i<n ' J 



(3.2) 



Proof. Set I n = E\(f>(X T ) - (j>(X£)\ 2 and for i = 1, ...,n, let 



h-i ■= 



E 



\Z r -Zl\ 2 dr. 



Using (fL3]) with £_ = 4>{X T ) and (|3~Tj) , we deduce 



Y 



-1 - + jf (^r - ) dW r =Y U - Yl + jf (/ (r, 9 r ) - / (fc, G^ 1 ) ) dr 

+ (. 9 (r,e r )-. 9 (t l ,6?'; 1 ))dB r . (3.3) 



DISCRETIZATION OF BDSDES 



7 



By construction, Y ti _ 1 — Y* is Tt i _ 1 = Ql i _ l measurable while for r e [ti-i,U), Z r — Z* is 
(Sr)-adapted. Hence, Y til — Y£_ is or thogonalto J t u _ (Z r - Z?)dW r . Therefore, 



I,-i =E 



Yu-l-Yu-i + / (Zr-Z?)dW r 



Since g(r, r ) (resp. g(U, Q^' )) is T r ( resp. J^J-measurable, the random variables Y ti — Y* and 
it* (g(r, X r , Y" r ) — g (tj,X^,l^)) dB r are orthogonal. Hence for every e > 0, using assumption 
[2J the L 2 -isometry of backward stochastic integrals, Schwarz's inequality and (|3.3p , we deduce 



/ (r,e r ) -/(*,-, GJ 1 ))^ 



1 + ^1 



^(l + Atie-^E^ -Y£| a + (Aii + 2e)E / / (r, 6 r ) - / [U, 0£ 



dr 



< 



+ (1 + A^e" 1 ) E / * g(r,e r )-g (t h QjA * dr 
1 + AUe- 1 + 2L f (At 2 + 2eAu) + 2L g (At, + At^e" 1 ) 



E \Y t . - Y f n \ 



L f (At, + 2e) +i 9 (1 + At.^ 1 
£/ (Atj + 2e) + a(l + At.e" 1 



E 



ti-i 



ir\ + \X r -Xl\ +2 



I^-FJ 2 ) 



dr 



Zf Z £ . 



dr. 



For |7r| < 1, At 2 < Ati] using Proposition 12 . 1 1 with p = 2 and Proposition 13. 11 we deduce 



E + \X r -Xl\ 2 + 2\Y r -r ti | 2 ) ^ < C|tt| 2 , 



for some constant C > 0. Hence for any 7 > 



h-x < 



1 + (e" 1 + 21/(1 + 2e) + 2L g (l + e -1 ) ) At 



E\Y ti -Y£ 



C 



L f (Ati + e) + L g (1 + Ai^" 1 ) 



+ (l+ 7 - 1 )[i / (At l + 2e) + a(l + At t e- 1 )]E / \Z r - Z t f dr 



+ (1 + 7) [L f (At, + 2e) + a (l + At^ 1 ) ] At, 



Zu - K 
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Lemma 13.41 yields for some positive constants C c , C £ , 7 and C ei7ji g, we have: 

i-U 

h-i < (1 + C e Ati)E\Y ti - Y*\ 2 + C e \ir\ 2 + C e . 7 E / \Z r - Z t f 

Jti-! 



+ k (1 + 7) (1 + 0) L f (AU + 2e) + a (l + A^e" 1 ) 



E 



\Z?-Z r \ 2 dr 



+ k (1 + 7) (1 + /T % (At, + 2e) + a (1 + Afce -1 ) 



E 



< (l + CeA^ElF^-F^I' + CeH^ + C^E / \Z T -Z t Xdr 

Jti-! 

+ k(1 + 7 )(1 + /3) [L / (At ! + 2e)+a(l + At i e- 1 )j / + \Z? - Z r \ 2 dr 



(3.4) 



Recall that a < - and let < 5 < 1 — k<x Then choose positive constants (3 and 7 small 
enough to ensure k (1 + 7) (1 + f3) a < 1 — Finally, let e > small enough to ensure that 
2k (1 + 7) (1 + fi) Lf€ < |. Then (|3.4[) implies the existence of C > such that for every i = 
1, ...,n - 1, 



+ ^E [ + \Z* - Z r \ 2 dr<(l + CAt t ) Ii + C\ir\ 2 + CE [ + \Z r - Z t f dr. (3.5) 
Using the discrete Gronwall lemma in |Z04| (Lemma 5.4 page 479), we deduce 



max <Ce CT E[I n 

0<i<n 



E 



\Z r — Zf.\ dr + \tt 



ti-i 



<Ce(|0(X t )-0(X?)| 2 + 2 



•Zr — Zt. 



\ 2 + \Z r - Z t f)dr+\ir\). 



Since is Lipschitz, Proposition 13.11 implies that E\4>(Xt) — </>(XJ)| 2 < C\tt\; thus Theorem 
implies 



max E\Y U -Y?\ < C\tt\. 

0<i<n 1 



(3.6) 



0<i<n-2 



Moreover, summing both sides of (|3.5[) over i from 1 to n — 1 and using Corollary 12.41 we obtain: 

+ [ \Z?-Z r \ 2 dr< J2 (l + CA^/i + CW 

6 Jt * l<i<n 

+ C ^ E / l^r-^tj 2 ^, 
l<i<n "^"i 

<C|tt| + ^ (l + CAti)^. 



Therefore, 



S 

I + -E \Z; - Z r \ 2 dr <C\ir\ + + C ^ A ** J ' 



Ki<n-1 



Since 6 < 1 — ko < 3, using (|3 .6[) we deduce 

- .. / - Z r | 2 dr <C|tt| +E / " \Z* - Z r \ 2 dr + C\ir\E I \Z? - Z r \ 2 dr. (3.7) 
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The equations (|1.3|) and (|3.ip imply 

r'n 



/ " (Z; - Z r )dW r =(Y t l -Y t J- {Y t l_ t - K tB _J 



+ (/(t ni X£,i£,0)-/(r,X r ,Y r) Z P ))dr 



+ /""(<? (t ni Xl,Y t l,0) - .g (r, X r , 3y, Z r )) dB r 

Jt n -! 

The L 2 -isometry, Schwarz's inequality, (|3.6p . Lemma l2~2l Propositions 12.11 and 13.11 



/■tn 

E / \Z? - Z r \ 2 dr <4E \Y t l - Y tn | 2 + 4E 



Y? , - K 



+ 4|tt|E / " 1/(^,^,3^,0) -/(r,X r ,Y;,Z r )| 2 dr 



tn-l 



+ 4E / " | 5 (i„, X£ , 0) - g (r, X r , Y r , Z r )\ 2 dr 
Jt n -i 

<C\ir\ + C\ir\ sup E (\X r - X T \ 2 + \X? n - X T \ 2 ) 

t n -i<r<t n 

+ C\tt\ sup E(\Y r -Y T \ 2 + \Y t l-Y T \ 2 + \Z r \ 2 ) 

t n -i<r<t„ 

<C\ir\. (3.8) 
For |-7r | small enough, we have C\ir\ < 5/6; thus (|3.7p and (|3.8|) conclude the proof. □ 

4. A NUMERICAL SCHEME 

In this section we propose a numerical scheme based on the results of the previous sections. 
First of all, given x G R d , s < t we set: 

X t (s, x):=x + (t-a)b (x) + a (x) (W t - W s ) . 

We clearly have X£ = X ti (ti-\, X^A for every i = l,...,n. Then, givenavector (xo, . . . ,Xi\Xi+i, 
R (i+i)d x R «- 4j get Xn+i = an d for i = 0, . . . , n - 1, let 

x* := (x , ...,Xi) , x l+ i := (x i+ i, ...,x n ) . 

Define by induction, the functions (resp. the random variables 

Uf, V? : R( l+1 ) d x n x R"^- 1 -> R) as follows: 

< (a; , x n ) := (f> (x n ) , (%o, x n ) ■= 0, 

and for i = 0, . . . , n — 1 let 

C7f (x'\ w, x l+2 ) :=< +1 (x'\ X u+1 (t h an) , x l+2 ) (4.1) 

+ /ftt+li -^t»+i (tii x i) ! M i+1 

(x i ,X tj+1 (ti,a;i),x i+2 ), 
( x \^* I+ i {U,Xi) ,x l+2 ) )A* i+ i, 
(x>,x I+2 ) :=ff^ i+ i,X 4i+1 ari) , uf +1 (x l ,Xj i+1 (ij, x t ) , x l+2 ) , 

< +1 (x%X ti+1 (^),x i+2 )), (4.2) 
< (x*; Xj+i) :=Ef7f (*\ w, xj+i) + x l+ iEV^ (x i , w, x i+1 ) , (4.3) 
wf (x i; ;x 4+1 ) :=— ^E (f/f (x>,x i+1 ) AW tj+1 ) 

+ -g+^E (C (x>,x i+1 ) AW ti+1 ) . (4.4) 
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Theorem 4.1. We have for all i = 0, ...,n 

Y£ =u*(X? o ,... ) Xl;AB ti+1) ...,AB t 



XZ;AB ti+1 ,...,AB tn 



(4.5) 
(4.6) 



Proof. We proceed by backward induction. For i = n, by definition Yj* = <f> (Xf^, so (|4. 5[) and 
(|4~6)) hold trivially. 

Suppose that the result is true for j = n, n — 1, • • ■ , i. The scheme described in (|3.1|) implies that 

= i£ + / (u, e^ 1 ) At, + .g (t n e^; 1 ) AB ti - jT z r w r . (4.7) 

To prove (|4. 5|) . we take the conditional expectation of (|4.7j) with respect to Tt i _ 1 = TQ t . V T$ T ; 
this yields 
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Efgft^e^; 1 ) AB U 



T 



U-i 
ti 



-E 



ZZdWr 



ti- 



T 



Using the fact that J t " Z£dW r is orthogonal to any J 7 t 4 _ I -measurable random variable, and the 
induction hypothesis we deduce: 



Y tt 



T 



E / [t h e^ )AU T ti _, +E[glU,,el 



Tt^ AS 



ti-i 



uf , . ., , X ti (tw , ) , AB U+1 ,...,AB tn 
ABi E(.g(ti, X ti (ti_i, J, < (X%, X£_ l5 X tj (t,_ 1; X£_ J, AB U+1 , ...,AB tn ), 



>?(xZ,...,XZ_ 1 ,X u {t i - 1 ,XZ_ 1 ),AB u+1 ,...,AB t 



1 are T ti _ 1 measurable while W*. — Wt.._, is 



Since all AB tl , I = 1, . . . , n and X£ k ,k = 0, . . . , 
independent of Tt i _ 1 ] we deduce (|4.5I) . 

To prove (|4.6[) . multiply (|4.T|) by AW ti = — Wt i _ 1 and take the conditional expectation with 
respect to Tt i _ 1 , this yields 



•Ft 



- E AWi 



ZZdW r 



Since £ T ti _ 1 and AW ti is independent of an d centered we deduce 



( Y U-i 



AW, 



T t , 



Furthermore, 



El AWti / ZldWr 



F 



=E 



--E 



0. 



ZZdr 



U-i 



ZZdr 



TuA =AuzZ:\ 



this completes the proof of (14.61) . 



□ 
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